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Abstract 

We discuss the weak-coupling BCS theory of a superconductor with the 
impurities, accounting for their anisotropic momentum-dependent potential. 
The impurity scattering process is considered in the t-matrix approximation 
and its influence on the superconducting critical temperature Tc is studied in 
the Born and unitary limit for a d^2_y2- and {d^2_y2 + s)-wave superconduc- 
tors. We observe a significant dependence of the pair-breaking strength on 
the symmetry of the scattering potential and classify the impurity potentials 
according to their ability to alter Tc- A good agreement with the experimental 
data for Zn doping and oxygen irradiation in the overdoped cuprates is found. 

PACS numbers: 74.20.-z, 74.62.-c, 74.62.Dh 
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I. INTRODUCTION 



Several experiments probing the effect of impurities or lattice defects on superconduc- 
tivity in the cuprates have been carried outEHil in order to get more insight into the sym- 
metry of the superconducting state. The most thoroughly studied defects are Zn and Ni 
substitutionsEFEl on the planar Cu sites and irradiation induced oxygen vacancieJii^ in the 
copper oxygen planes. Yet the results are difficult to explain within a standard Abrikosov- 
Gorkov type theory of impurity scattering!! for the scenario of the ci-wave superconductivity, 
which is predicted to be extremely suppressed by the impurities.ii"0 This issue was critically 
examined by Radtke et al.@ who considered isotropic nonmagnetic impurity scattering in the 
Born approximation and obtained the critical temperature in both weak- and strong-coupling 
approach close to the Abrikosov-Gorkov scaling function. A comparison to the electron ir- 
radiation data@ in Y Ba2Cu^0-j^s (F — 123) showed that the theoretically predicted was 
about twice as much reduced by the impurities than observed. This inconsistency can be 
settled down within the weak impurity scattering model if the impurity scattering rate a fac- 
tor of 3 less than the one deduced from the transport measurements is assumed^S which is 
equivalent to an introduction of two separate relaxation time scales - one defining the scatter- 
ing time in pair-breaking processes and the other representing the transport scattering time. 
Such a distinction occurs naturally as a consequence of an impurity momentum-dependent 
scattering probability.il The issue of a possible anisotropy in the impurity scattering poten- 
tial was suggested in a discussion of the irradiation data by Giapintzakis et al.lii where the 
authors evoked a model by Millis et al.il A more general formulation of the problem was 
given in our previous paperji^ where the effective correlation between two impurity vertex 
functions was assumed in the form |u;oP + l^^iP/ (k) / (k'), with \wq\ and \wi\ representing 
isotropic and anisotropic impurity scattering amplitude respectively and / (k) determining 
the symmetry of the impurity potential.0 Our analysis showed that the symmetry of the 
anisotropic potential is an important factor and a significant reduction in the pair-breaking 
strength appears for large values of (e/)^ = [/^^ dStn (k) e (k) / (k)]^, where jpg dSk denotes 
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integration over the Fermi surface (FS) and n (k) is the normahzed {^ps dSkU (k) = 1) angle- 
resolved FS density of states. Particularly low T^. suppression is predicted for (e/)^ = 1 that 
is for the anisotropy of the scattering potential in phase with the order parameter, which 
in the case of the dx-2-y2-wme superconductor corresponds to / (k) ~ — ky. Within this 
approach we were also able to understand quantitatively the irradiation data@ in F — 123. 

Through the postulation of the analytic form of the square value of the impurity potential 
this approach has been designed for the second-order Born approximation and cannot be 
efficiently extended to include multiple impurity scattering processes. Such a generalization 
is important from the theoretical and experimental points of view, especially in the light 
of the recent experiments suggesting a possible strong (close to unitary) scattering by Zn 
atoms in F — 123 and La2-xSrxCu04 {La — 214) compounds.! For that purpose a model 
based on the assumption of the impurity potential and not its square value is needed. 

In this paper we study in the t-matrix approximation the pair-breaking effect of the 
nonmagnetic impurities with the anisotropic momentum-dependent factorizable potential. 
The influence on the superconducting transition temperature is analyzed quantitatively in 
the Born and unitary scattering limits. In particular, we flnd that the scattering from 
impurity potential in phase with the order parameter leads to a stronger suppression of the 
critical temperature than from the impurity potential orthogonal to the superconducting 
state. On the other hand, the superconducting state appears very robust to the potential 
scattering with its maxima in the region of the nodes of the order parameter. Finally, we 
flnd that the Zn and electron irradiation suppression data in F — 123 and La — 214 are in 
the range predicted by our model. 

This work is organized as follows. In Sec. II, we introduce the anisotropic momentum- 
dependent impurity potential. In Sec. Ill, we derive the expressions necessary for the 
analysis of the scattering process of arbitrary strength. In Sec. IV, we obtain the self- 
energies due to impurity scattering at the superconducting-normal state phase transition 
which allows the evaluation of the critical temperature. Then we examine analytically Born 
and unitary scattering limits. In Sec. V, we calculate numerically the critical temperature 
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for d- and {d + s)-wave superconducting states in the presence of impurities accounting for 
their anisotropic potential of p-, d-, /-, and g-wave symmetry in Born and unitary limit. 
In Sec. VI, we compare the results with data on Zn substituted and irradiated overdoped 
La — 214 and optimally doped F — 123 samples. In Sec. VII, we present our conclusions. 
Except for comparing to the experimental data we assume /i = /c^ = 1 in the calculations. 



As the impurity scattering strength is rather impossible to be determined from first 
principles, the information about the scattering process is usually deduced from a compar- 
ison of the experimental data with the theoretical models. These models assume a certain 



data. We proceed in the same way by analyzing the impurity scattering potential of the 
factorizable form 



The above interaction consists of two channels - the isotropic scattering channel with the 
scattering amplitude Vi, and the anisotropic one determined by the scattering strength Va 
and momentum-dependent function / (k). We assume that the / (k) average value over the 
Fermi surface (/) = jpgdSkn (k) / (k) = 0. Therefore, / (k) is orthogonal to the isotropic 
s-wave term in a sense of a scalar product defined as the FS integral and in consequence 
a symmetry other than the identity is introduced into the impurity potential. An addi- 
tional normalization (/^) = 1 gives Va the meaning of the scattering strength magnitude in 
the anisotropic channel. The chosen potential depends on the absolute orientations in the 
crystal of the incoming and outgoing (scattered) particles momenta, not only on the angle 
between k and k'. The usually assumed isotropic condition^ of a spherical (cylindrical 
in two dimensions) constant energy surface and the scattering probability dependent only 
on the angle of deflexion are broken here.S This leads to a momentum-dependence of the 



II. IMPURITY SCATTERING POTENTIAL 





(1) 
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time between scattering events determined by the imaginary part of the self-energy and a 
momentum-dependent relaxation time in the Boltzmann equation.H'il This second quantity 
is worth mentioning as it provides the information about the normal state transport proper- 
ties and can be used as an additional physical assessment of the phenomenological model.ii 
A constraint on the potential (1) which follows immediately from the analysis of the normal 
state properties is a nonzero value of the scattering amplitude in the isotropic channel fj. 
A lack of the s-wave scattering may lead in some cases to an infinite value of some elements 
in the dc conductivity tensor,ii that is to nonphysical transport properties. 

The anisotropic potential from Eq. (|l|) distinguishes a certain set of coordinates which we 
think should coincide with the main directions in the crystal. This assumption seems rather 
plausible because the impurity potential is determined by the dopant atom itself as well 
as its substitution site in the crystal. The potential produced by a given sort of impurity 
(defect) may be considered unique since the impurities tend to selectively substitute at 
characteristic sites in the crystal. The Zn and Ni atoms occupy the in-plane Cu sites0^ii 
and the electron irradiation displaces the oxygen atoms from the Cu02 planes.liiHHi Therefore 
on a short length scale of the order of magnitude of the lattice constant a given impurity 
dopant produces the same potential of the same orientation throughout the crystal. On the 
other hand, on a large length scale the impurity distribution in the system is random and 
the Abrikosov-Gorkov's method^! of averaging the Green's functions is applicable. 

We have not made any assumption about the electron energy band yet and we use 
a general formalism deriving the equations valid for an arbitrary Fermi surface. For the 
computational simplicity, however, the numerical results are obtained for a cylindrical FS. 
This approximation allows to study the effect of the anisotropy of the impurity potential 
alone. The anisotropy of the Fermi surface may equally enhance or suppress the pair- 
breaking effect of the impurities. The detailed quantitative calculations are needed to answer 
this question which is beyond the scope of the present paper. 

For the numerical calculations we take the function / (k) in the scattering potential (|Tp 
proportional to the harmonic functions which in a polar angle notation read sin{l(j)) and 



cos{l(f)), where / is an integer number. Therefore we study the effect of the basis elements 
in a space of the functions determined by a two-dimensional momentum vector. 



III. T-MATRIX APPROXIMATION FOR THE SELF-ENERGY 

We study the effect of potential scattering by spinless, noninteracting impurities on the 
single-particle propagator for superconducting electrons 

G (k, Lu) = iujfoao - ^kTscro - A (k) iT2a2 - S (k, u) (2) 

Here C,k is the quasiparticle energy, u = 7rr(2m +1), where T is the temperature and m 
is an integer, fj, crj {j = 1,2,3) are the Pauli matrices and tq, ctq are the unit matrices in 
particle-hole (Nambu) and spin space respectively. The order parameter A (k) is defined as 

A (k) = Ae (k) (3) 

where e (k) is a momentum-dependent real function which may belong to a one-dimensional 
(ID) irreducible representation of the crystal point group or may be given by a linear com- 
bination of the basis functions of different ID representations. We normalize e (k) by taking 
its average value over the Fermi surface (e^) = 1. This normalization gives A the meaning 
of the magnitude of the order parameter. The self-energy S (k, uj) and consequently the 
Green's function (Eq. (^) have been obtained by applying Abrikosov-Gorkov's techniqueil 
of averaging over the coordinates of the impurities and depend on only one momentum 
vector k. In this approximation 

±{k,uj) =nf{k,k,uj) (4) 

where n is the impurity concentration and T obeys the Lippmann-Schwinger 
equationii'ii'&^i 

f (k, k', u})=v (k, k') + J2v (k, k") G (k", u) f (k", k', u) (5) 

k" 

Since the scattering potential for a single electron 



v{k,k')=v{k,k')f,ao (6) 

is momentum-dependent (Eq. (|l])), the vertex part T (k, k', uj) is a function of two momenta 
k and k'. Therefore it should be evaluated from Eq. using the explicit form of v (k, k') 
first, and then the self-energy can be obtained according to Eq. by taking k' = k. We 
proceed to a solution by defining 

gi{uj)=Y.r{k)G{k,uj) , 2 = 0,1,2 (7) 

k 

and expanding all matrix quantities as 

S (k, uj) = So (k, Lo) foo-Q + Si (k, uj) fio-2 + S2 (k, u) + S3 (k, u) t^gq (8) 

cji {uj) = Qio (uj) fo&o + Qii (uj) ficr2 + gi2 (u;) r2a2 + ga (u;) Tso-q (9) 

The expression for the one-particle Green's function is then 

A X _ i^To&o + 4^30-0 + A (k) ifso-a + A' (k) fiaa , . 

U[K,UJ — ^2 ~ ~ \^^) 

u^ + ^k - A2 (k) + A'2 (k) 
with (D = LU+iTiQ (k, u), = C,k+^3 (k, cj), A (k) = A (k)— ■iS2 (k, u) and A' (k) = Si (k, uj). 
The formal solutions for the self-energies S^ (k, uj) (j = 0, 1, 2, 3) obtained through Eqs. (H) 
and d^) read 

Sj (k, uj) = Sj {vi, Va, go, 91,92) u {vi, Va, 90, 91,92) 

+ Sj {Va, Vi, 92, gi,go) U {Va, Vi, g2, (ji, 9o) P (k) (H) 

+ [t-i_j {vi, Va, 90, gi, g2) + h~j {va, Vi, g2, gi, go)] f (k) 
Note a permutation: Vi ^ Va, §0 ^ §2, in the arguments of the anisotropic parts (propor- 
tional to / (k) and (k)) of Sj {k,uj) in Eq. (|ll|). The functions sj, tj and u are given by 
a series of equations 

So {vi, Va, go, 91, 92) = floCs + aiC2 + a2Ci + asco 
siivi, Va,go, 91,92) = 0.0C2 + aiOi + ia2Co- ia^ci 

(12) 

52 {vi, Va, go, 91,92) = floci - iaico + a2C3 + m3C2 

53 {vi, Va, go, 91,92) = floco + iaici - ia2C2 + 0303 



tj (vi, Va, go, gi, h) = ViUCj (cg + cl + cl- cfj ^ (13) 

u (fj, Va, go, gi, h) = ViV^^nd'^ (cq + cl + cl- Cg) (14) 

with d = — (710 + +5'i2 + fi'i3 the coefficients aj, Cj (j = 0, 1, 2, 3) determined by twelve 
gik elements (Eqs. (|^) and (|)) given by the integrals of the products of the Green's function 
G (k, a;) and appropriate powers of the impurity potential anisotropy function / (k). These 
coefficients are introduced in order to shorten and simplify the notation of the self-energy 
functions. We define them in the sequential formulas with a use of additional bj parameters 

OO = fl'13 (1 - Vag23) " Va (fi'l2fl'22 + fi'llfi'21 " fi'lo5'2o) 
ai = igi2 (1 - Vag23) + Va («5'l35'22 " 5'll5'20 + 5'l0fi'2l) 

0-2 = -igu (1 - fafi'23) - Va {igizg2i + gi2g2o - giog22) 

as = -giO (1 - Vag23) - Va (fi'l3fi'20 " «fi'l2fl'21 + «5'll5'22) 



(15) 



(16) 



(17) 



60 = (1 - Vigos) oo + Vi {igo2ai - ^01^2 - 5^00^3) 

61 = (1 - Vigos) ai + Vi {igo2ao + igooa2 + goi^s) 

62 = (1 - Vigos) a2 + Vi {go2a3 - Woi^o - igoo^i) 

63 = (1 - Vigos) as - Vi (5^0202 + goiai + gooao) 

Co = v^'^d-% - Vigis 
ci = iv~'^d'% - Vigi2 
C2 = -iv~^d-% - Vigil 
C3 = -v-^d~% + v^gio 
The self-energies can be evaluated with a simultaneous solution of the gap equation 

A (k) ia2 = -Tj2J2 ^^k,k'tr [{n + G (k', to)] (18) 

where Vk,k' = — Vqc (k) e (k'), Vq > 0, is the pair potential. The above lengthy expressions 
have been derived without any additional constraints and are fundamental to the considera- 
tions within the model. This generally complicated problem simplifies with the assumption 

8 



of particle-hole symmetry of the excitation spectrum. It has been shown by Hirschfeld et 
al.Ell that in this case the 7^3 component of the integrated Green's function ((703) may be 
neglected in the presence of s-wave scatterers. Inclusion of higher-order angular momentum 
waves in the scattering potential makes this analysis considerably more difficult. In this pa- 
per, however, we calculate the self-energies at the phase transition and for that purpose we 
need to show only the consistency of the assumption S3 = at the critical temperature Tc. 
In the following we assume a particle-hole symmetry of the energy spectrum, take (^43 = 
(i = 0, 1, 2) and check if this condition leads to a vanishing self-energy E3. 



IV. SELF-ENERGY AT PHASE TRANSITION 



We consider the effect of anisotropic impurity scattering on the critical temperature. At 
the superconducting-normal state phase transition the gap equation (O) transforms into 



_ _ e (k) + 1 (El (k, uj) /A)^^o + - (S2 (k, uj) /A)^^o 

UJ k 



{uj + iSo (k, ^)a=o)^ + i^k + S3 (k, ^)a=o)^ 



(19) 



In order to find Tc the self-energies in A — limit are to be obtained. For the sake of 
convenience we introduce new parameters q = l/(7riVofj) and Ca = l/i^nNoVa) describing 
the scattering strength in the isotropic (cj) and anisotropic (c^) channel respectively. A^o 
represents the overall single-spin density of states at the Fermi level. A new measure of 
impurity concentration F = n/(7rA^o) is also used. Taking into account that (fi'oo)A=o ~ 
—mNQsgnioj) and (5'2o)a=o ~ (fl'oo)A=o (because of normalization {p) = 1) we obtain 

1 1 



A=0 



-iT 



ci + 1 cl + l 



f (k) 



sgn [UJ) 



(20) 



S3 (k,tc;) 



+ 



-f (k) 



(21) 



^ [cl + l ■ cl + r 

Calculation of (Si (k, cj) /A)^^q and (S2 (k, cu) /A)^^q quantities is more tedious and re- 
quires a solution of two sets of two linear equations which determine four unknowns 
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(5'ji/A)^^Q, (fi'j2/A)^^p (j = 0,2) (note that gn=gi2 = because of (/) = 0). This procedure 
leads to (Si {k,uj) /A)^^q = and 



r 



where 



and 



^ (^7o2/A)^=o + (922/ A) f (k) 



sgn{uj) (22) 



TTiVo 



c? + i-r(- 



(^722/A)^^0 = -t [cl + 1) 



(23) 



X 



ci + i-r 



c.' + 1 - r 



(24) 



X 



c +1-T 



1 



with ljq = u + iTiQ (k, ti^)^=o- following subsections we consider the above self-energies 

in the Born and unitary scattering limits. 



A. Born scattering 



When both isotropic and anisotropic impurity scattering channels are in the Born scat- 
tering regime i.e. q 1 {viNo <^ 1), Ca ^ 1 {vaNo <^ 1) only the lowest order terms in the 
impurity potential play role in the self-energies. Keeping up to the square terms in Vi [va) 



we obtain from Eqs. (EOf) and (BTI) 



iirnNo + vj (k) sgn (uj) 



(25) 



S3 (k, uj)^^Q = n Vi + Vaf (k) 



(26) 
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Though S3 is nonzero, it may be absorbed into the chemical potential and its effect vanishes 
in the Born scattering.il For (S2 (k, lu) /A)^^q we obtain 

(E2 (k, Lo) I A)^^Q = mnNQ 

It is noteworthy that apart from the average value of the order parameter symmetry 
function, (e) , a term reflecting the overlap between e (k) and p (k) influences the self- 
energy (S2(k, cj) /A)^^Q. The limit of fa = gives the standard s-wave impurity scat- 
tering in the Abrikosov-Gorkov approximationil with So(k, C(j)^^q = —i7mNQvfsgn{uj), 
(S2 (k, cj) /A)^^Q = innNovf {e) sgn (lj) /u, and the critical temperature determined by 
In (TJTJ = {{ef - l) (1/2 + nN.vf/ (2T,)) - ^ (1/2)]. 



{e) + vUef)f{k) 



(27) 



B. Unitary scattering 



The limit of the resonant impurity scattering in both isotropic and anisotropic channels 



i.e. Cj — 0, Ca — > in Eqs. (EOf) and (pl|) leads to 



So(k,a;)^^o = -zr 1 + f (k) 



sgn [uj) 



(28) 



and S3 (k, Ci;)^^Q = which is consistent with the assumption of particle-hole symmetry of 



the excitation spectrum. The unitarity limit in Eqs. (|2^ ) - (|2^) gives 



(E,(k..,/AU.„^>r{[i-r(-l)]"(^) 



+ 



X 



i-r(± 



1 -r 



(29) 



\UJq/ / \UJ0 \UJ0 \UJ0 



sgn [uj) 
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The impurity potential may lead to the strong scattering in one channel and the weak 
scattering in another. If the strong scattering takes place in the isotropic channel and the 
scattering in the anisotropic one is in the Born limit, then we deal with a case of the isotropic 
unitary scattering.iiS&^i The opposite case with the strong scattering in the anisotropic 
channel and the weak scattering in the isotropic one is equivalent to a nonphysical situation 
of unitary scattering in the anisotropic channel alone as discussed in section II. 



V. CRITICAL TEMPERATURE 

The critical temperature in the weak-coupling BCS approximation is determined by the 
equation 



T 

, CO , 



sgn {u) 



1 _ e(k) + l(Si(k,a;)/A)^^o + -(S2(k,a;)/A)^^o 
E e (k) ' ' 



TTiVo V + iT.^ (k, u)^^^Y + + S3 (k, u)^^^y 

(30) 



where because of the momentum-dependent self-energies the summations over the quasimo- 
mentum vector k restricted to the Fermi surface and Matsubara frequency uj need to be 
performed numerically. We do the calculations in two different scattering regimes - Born 
and unitary for the rf-wave order parameter. Finally the results for the {d + s)-wave super- 
conductor are discussed briefly. In order to proceed further one has to choose a function 
describing the anisotropy of the impurity potential. We study the functions given by the 
subsequent harmonics as they represent an ortho normal and complete set in the interval [0, 
27r] and can be used for the Fourier expansion of any regular function. Harmonics up to the 
4th order are analyzed and the results are extended to higher-order harmonic functions. It 
is important to realize that the impurity potential of p-wave anisotropy shows a particular 
property. Given by sincf) or coscj) function it represents a small-angle scattering which is 
relatively unimportant in contributing to the resistivity of the normal state. Taken even to- 
gether with the s-wave scattering channel it results within the linear response approximation 
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to the distribution function in a not well defined integrals determining the dc conductivity 
for some electric field orientations.il In other words, in the Boltzmann equation based anal- 
ysis of the transport properties the p-wave scattering amplitude couples to the "1 — cos(6')" 
term in the collision integral. Therefore, in a fit to the real systems the p-wave scattering 
can be considered only as one of the components in the anisotropic channel of the impurity 
potential coexisting with some other higher-order harmonics. In order to emphasize this 
different feature of the p-wave potential among the other basis functions we present the 
momentum-dependent scattering with various amounts of the p-wave potential in a separate 
figure. 

The anisotropic impurity scattering is compared to the isotropic one of the same scatter- 
ing strength. It means that we discuss the effect of substituting a part of a given isotropic 
impurity potential Vq with an anisotropic term. The amplitude of the replaced isotropic 
potential, fo, and the amplitudes of the isotropic and anisotropic scattering channels in the 
studied potential are related through the formulas Vi = avo, and = (1 — a) Vq, where the 
coefficient a (0 < a < 1) defines the partition of the potential. 

A. Born scattering 

The impurity scattering is analyzed as a function of the pair-breaking parameter T' = 
ttuNoVq (fo = Vi + Va) which includes the overall scattering strength that is, takes both 
isotropic and anisotropic channels into account. It is convenient to use the ratio Va/vi = 
(1 — a) /a to define a particular potential. 

1. s-wave superconductor 

We note, before discussing the unconventional superconductivity, that for the isotropic 
s-wave superconductor, given by e (k) = 1, the Tc equation ( pO| ) along with Eqs. ( pSj ) and 
(P7|) do not lead to a change of the critical temperature in agreement with the Anderson's 
theorem.B 
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2. d-wave superconductor 



The critical temperature for a dj.2_j^2-wave superconductor that is for e (k) ~ (k^ — ky) 
is presented for a large range of Va/vi ratio values in Fig. 1. The major common feature of 
these diagrams is a lower Tc suppression by the anisotropic impurity scattering compared to 
the isotropic (s-wave) one. According to their effect on superconductivity the anisotropic 
potentials can be classified in three groups defined by the function / (k) : cos2(j), sin2(f), 
and higher-order harmonics. The first of them, determined by / (k) ~ cos2(j), leads to the 
strongest Tc suppression. In this case the directions of the maximum impurity scattering 
correspond to the maxima of the order parameter and most of the pair-breaking process 
takes place in this region. Therefore, the suppression of superconductivity is particularly 
strong. However, for the level of anisotropy up to (1 — a) 0.5 {va/vi 1) the depairing 
effect of impurities is reduced with increasing contribution of cos2(f) scattering compared to 
the isotropic scattering. Further increase of the anisotropic part in the impurity potential 
up to a level (1 — a) 0.6 {va/vi ^ 1.5) almost does not change the pair-breaking effect. 
When the contribution of the anisotropy exceeds this level the impurity effect on Tc is 
enhanced and finally for (1 — a) of the order of 0.83 {va/vi = 5) becomes comparable to 
the one of the s-wave scatterers as it is shown in Fig. le. The out of phase scattering 
takes place for / (k) ~ sm20. This function is a 45° rotation in the xy plane of the order 
parameter function e (k) ~ cos24>, that is, the impurity potential maxima correspond to 
the superconducting gap nodes and vice versa. In this way the impurity pair-breaking 
effect is minimized. The suppression of the critical temperature is reduced by an increasing 
amount of the anisotropic scattering in the impurity potential. For large Va/vi, however, 
the shape of the suppression lines changes slightly and a small enhanced suppression can 
be observed for some impurity concentration (Fig. le). In general, the anisotropic impurity 
potential given by a function orthogonal to the order parameter is less pair-breaking than 
/ (k) ~ e (k) . The Tc suppression by the 3rd and 4th order harmonics is less than the one 
of cos2(f) but it exceeds that of sin2(f). It is also decreasing for increasing level of anisotropy 
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in the impurity potential up to (1 — a) ~ 0.5, then the pair-breaking effect of impurities 
practically saturates. The curvature of the graphs changes at some points, however, and the 
normalized critical pair-breaking parameter, V / (27rTc(,), at which Tc ~ may increase in 
some cases. It is worth observing that the critical temperature is almost equally suppressed 
by the impurity potentials given by the 3rd and 4th order harmonics for Va values up to 
Va/vi = 1. This is particularly true for small anisotropic scattering levels. When the 
contribution of the anisotropic channel in the impurity potential is larger, the differences 
are more pronounced. Nevertheless, even for the amplitudes ratio as high as Va/vi — 2 
( 1 — a ~ 0.67 ) the curves differ only near a zero value of the critical temperature and the 
normalized critical pair- breaking parameters, V / (27rTc(,), of these four harmonics are within 
an interval of the order of magnitude of 10~^. For the sake of transparency we present cos3(f> 
scattering effect only in Fig. Id. Though the curves of different harmonics overlap for the 
most of the critical temperature range for the amplitude in the anisotropic scattering channel 
five times as much as the one in the isotropic channel, Va/vi — 5, they split distinctly at a 
low temperature of about 0.2 T^, (Fig. le). A study of the higher order harmonics up to the 
10th order for the anisotropic scattering strength equal to the isotropic one {va/vi — 1) leads 
to almost the same result as the 3rd and 4th order harmonics. This approximately universal 
behavior for harmonics from the 3rd to the 10th order is shown in Fig. 2. Concluding 
we may say that to a good accuracy the anisotropy of the impurity potential given by the 
harmonics of the order higher than two yields an approximately universal Tc suppression. 
It is also worth mentioning that except for the isotropic scattering the critical temperature 
goes to a zero value asymptotically. The asymptotic tails start at very low temperatures 
of the order of magnitude of 10~^ Teg (not seen at the figure scale) where the fluctuation 
effects become important and may destroy superconductivity. 

Remembering that this potential cannot be considered on its own without any higher- 
order admixture, we show the pair-breaking effect of the p-wave scattering in the anisotropic 
channel for different Va/vi values in Fig. 3. Compared to the other harmonics (Figs. 1-2) 
it yields the lowest Tc suppression. This fact can be explained by an effectively small- angle 
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scattering of the j9-wave potential. If it is given by cos0 function then mostly the quasi- 
particles with their momenta parallel to the x-axis are affected, and for the sincj) function 
representing impurity potential the electrons moving along the y-axis are being scattered. 
Thus we deal with a weak practically one-dimensional scattering in a two-dimensional space. 
The pair-breaking effect of both basis functions cos0 and sinip is the same within the ac- 
curacy of the numerical calculations. Worth observing is also a feature of a reduced 
suppression with an increasing amount of p-wave scattering in the impurity potential. How- 
ever, at large amplitudes in the anisotropic channel {va/vi ~ 5) the initial suppression of 
the critical temperature can be enhanced. Characteristic is also a clear asymptotic decrease 
of the critical temperature to its zero value. 

Any impurity potential given by Eq. (|l]) can be represented as a combination of the 
potentials based on single harmonics. Therefore, we expect that the pair-breaking effect will 
be given by an appropriate superposition of the effects discussed above. 

3. (d+s)-wave superconductor 

As an example of a {d^-i-y-i +s)-wave superconductor we consider the one with an s-wave 
admixture of 10%. This is the order of magnitude of the s-wave level which cannot be 
ruled out by the ARPES measurements in the cuprates.SS The normalized to unity order 

/ 2 \ 1 /2 

parameter is given then by e (k) = (cos20 + s) / ({cos2(j) + s) ) , where s = (0.005/0.99) ' 
and its FS average, (e) = 0.1, corresponds to 10% of the s-wave fraction in the {d^2_y2 + s)- 
wave superconductor. The results for Va/vi = 1 are shown in Fig. 4. Even a small amount 
of the s-wave component, which is not destroyed by the potential scattering, results in a 
robustness of Tc compared to the pure d-wave superconductor (Fig. Ic) and its asymptotic 
reduction. A change in the symmetry of the order parameter also causes a split in the 
suppression of the critical temperature by the j9-wave anisotropic channel scattering. 
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B. Unitary scattering 



The relation determining the partition of the scattering amphtudes between the isotropic 
and anisotropic channels, fj + Vaf (k) / (k') = v^a + Vq {1 — a) f (k) / (k'), holds as long as 
these amplitudes are finite. In the case of the resonant scattering we make an approxima- 
tion of Vi and Va diverging to infinity. Therefore we cannot control the relative scattering 
strengths in both channels and these processes become independent. Each part of the im- 
purity potential enters the self-energy functions (Eqs. through variables Cj and Ca- 
It is important to note that the self-energies depend on Cj and parameters separately 
so they are functions of q and Ca and not of any combination of them (like Ca/ci for in- 
stance). In the unitarity limit q — 0, Ca — > and the contributions from the isotropic 
and anisotropic channels are in fact equal. Similarly to the Born scattering limit, we want 
to discuss the effect of an anisotropic impurity potential that replaces the isotropic one. 
In order to do an appropriate comparison we need to refer to the isotropic scattering in 
two channels or equivalently to a regular isotropic scattering in one channel with a doubled 
impurity concentration 2n. Therefore, to compare the anisotropic unitary scattering with a 
corresponding isotropic one, the impurity concentration for the anisotropic scattering model 
must be half of the concentration of the s-wave impurities. We take the impurity concen- 
tration equal to n/2 for the two-channel anisotropic potential and the results are plotted as 
a function of F' = 2r = n/ (ttNo), where n is the real impurity concentration the same as in 
the Abrikosov-Gorkov scaling function for the s-wave scattering. That is, the two-channel 
scattering potential is averaged over two impurities so there is only one scattering channel 
per impurity present. This procedure introduces the anisotropic scattering potential of the 
form given by Eq. (P and the same scattering strength as the isotropic one in the unitarity 
limit. 
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1. s-wave superconductor 



The lack of influence of the potential anisotropic impurity scattering in the unitary 
limit on the critical temperature can be shown rigorously for a small concentration of the 
defects n/(7rA^o) < 1, when the self-energies So (k)^^^ = -iV [1 + f (k)] sgn (uj) (Eq. (|2|)) 
and (S2(k)/A)^^Q = — Eq (k)^^^ /tu (Eq. (^^). For a larger impurity concentration a 
numerical analysis confirms this result with a very good accuracy. 

2. d-wave superconductor 

The effect of the anisotropic unitary scattering on a dx2-y2-waMe superconductor is shown 
in Fig. 5. The same general rule as for the Born scattering holds here. The weakest pair- 
breaking effect is caused by the impurity potential involving sin2(j) and most suppression is 
seen for the potential containing cos20, that is the anisotropic scattering channel in phase 
with the order parameter. All the other harmonics orthogonal to cos20 yield a moderate 
suppression of the critical temperature comparable to the isotropic scattering. They also 
form an almost universal suppression curve. A study of the Tc reduction by the harmonics 
of the order up to ten shows the same approximate universality as that in Fig. 2, with 
the differences in the normalized critical pair-breaking parameter, F'/ {2ttTcq), of the order 
of 10"'^. The curves of the 3rd and 4th order harmonics overlap in Fig. 5. Particularly 
interesting feature of the resonant scattering is a strong pair-breaking effect of the cos2(f) 
anisotropic potential. It destroys superconductivity even faster than the isotropic impurities. 
Another characteristic fact is practically the same critical temperature dependence on the 
impurity concentration (i.e. not distinguishable in Fig. 5) for the p-wave scattering in the 
anisotropic channel and sm20 scattering potential. Similarly to the weak scattering limit, 
except for the s-wave scattering the suppression of superconductivity at low Tc is asymptotic, 
however not seen at the scale of Fig. 5. 
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3. (d+s)-wave superconductor 

A {dx2_y2 + s)-wave superconductor is more robust against the impurity scattering due 
to a nonzero s-wave component. Since cos20 function is no longer in phase with the order 
parameter its pair-breaking effect is lowered and becomes even less than the isotropic one for 
low critical temperatures. The results for the same level of the s-wave part in a {d^2_y2 + s)- 
wave order parameter as the one discussed in the Born limit are shown in Fig. 6. It is worth 
mentioning that the p-wave scattering leads again to the same Tc suppression as the sin2(j) 
anisotropic potential. 



VI. COMPARISON TO EXPERIMENT 



In the overdoped samples Tc/Tc^ data points plotted vs. impurity concentration form 
a universal curve independent of the critical temperature in the absence of impurities.ll2E 
Such a universal scaling behavior is characteristic of the impurity limited superconductivity 
provided the pair-breaking parameter F'/ (27rTc„) does not change with Tc^. This requirement 
is equivalent to a constraint NqTc^ = constant in the unitarity limit and may be obeyed 
in the overdoped systems where the critical temperature decreases with increasing hole 
concentration.0 There are also indications that a large residual resistivity due to Zn atoms 
in the cuprates corresponds to an impurity potential scattering in the unitary limit.! Thus, 
we consider the case of the unitary scattering focusing on the effect of Zn dopant. Working 
in the Born scattering limit requires an estimation of the impurity scattering potential which 
can be obtained from the residual resistivity in the normal state.@@iilii@ In that analysis 
a clear distinction between the impurity scattering life time and the transport relaxation 
time, which may differ if the impurity potential is anisotropic,!! is needed. In order to 
compare our unitary scattering results with the experiment we have to convert the pair- 
breaking parameter F'/ (2ttTcq) into the impurity concentration. As F' = n/ (ttNo) only 
the values of the density of states on the Fermi surface and the critical temperature in the 
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absence of impurities are needed for that purpose. The density of states at the Fermi level is 
estimated from the measurements of the specific heat jump at the phase transition AC.0E1 
We employ the BCS weak-coupling relation AC/7TC ~ 1.43 to obtain the normal-state 
Sommerfeld constant 7 which gives the density of states through 7 = 27r^A;^A'o/3. It is 
important to note that the strong-coupling corrections^'0 as well as the interaction with 
impuritie^ilii may change this relation significantly and in consequence alter the overall 
density of states. It may result in a wide range of Tc solutions.ii Because of a difficulty in 
the separation of lattice and electron contributions to the specific heat the thermodynamic 
experiments provide AC values with the accuracy depending on the quality of a sample. For 
high purity, fully oxygenated F— 123 compound grown in BaZrO^ the mean-field component 
of the electronic specific heat jump is estimated as 56 ± 2 {mJ / K'^mole). ESil In La-214 
system this quantity is in the range of 14 ± 5 {mJ / K'^mole)^^ There is no specific heat 
jump observed at the phase transition in Bi— 2122 compouncS thus the present method of 
obtaining the density of states cannot be applied in this case. We discuss the experimental 
results in this system together with F— 123 and La— 214 compounds more extensively in Ref. 
53. The density of states calculated from the above AC values are used for the evaluation 
of the impurity concentration, which is represented as the per cent number of impurities 
(defects) per planar Cu site. We do the calculations with a fixed density of states as we 
do not have any quantitative data showing its change with doping or disorder. Therefore, 
the results are not really universal and depend on the critical temperature in the absence 
of impurities T^q. We group them according to the values of the critical temperature of a 
pure system. In the optimally doped F— 123 compound T^^ varies in a very narrow range of 
values and for the theoretical calculation we take Tc^ = 91K. The phase diagrams for the 
d-wave superconductor together with the experimental data of Zn doped samples^Ull are 
shown in Fig. 7. The scattering effect of a given potential corresponds to the area between 
two curves of the same style, for instance the solid lines stand for the sm2(/) impurity 
potential. This broadened range of values stem from the uncertainty in the estimation 
of the FS density of states, that is, from the accuracy of the AC values. We note, that 
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all the experimental points fall in the region of the theoretically predicted Tc suppression 
by the unitary impurity scattering. Interestingly, the pair-breaking effect of Zn atoms in 
certain samples corresponds to the resonant scattering in different anisotropic channel of 
the impurity potential. We can clearly distinguish data that can be approximated by the 
scattering potentials given by cos20 and sin2(f) functions. The rest of the experimental 
points lie in the range of isotropic scattering or anisotropic scattering with higher order 
harmonics. Although the critical temperatures of the analyzed samples are very close, we 
do not see a universal suppression dependence as in Bi— 2122 and La— 214 compounds.0lil 
This feature may suggest a possible sample-dependence of the effective impurity scattering 
caused for instance by the differences in a sample preparation or Zn substitution processes. 
Even small changes in the hole concentration alter the electronic density of states at the 
Fermi level and result in a modified impurity scattering rate. 

According to Teg values, we gather the experimental data for Zn doped La — 214 in two 
groups corresponding to equal to 30K (Fig. 8a)0 and 367^ (Fig. 8b).0'i For the sample 
of Teg ~ 30K the experimental suppression of the critical temperature is in the range of 
anisotropic scattering determined by sm20 (Fig. 8a). The pair-breaking effect of Zn atoms 
in the samples characterized by Tc^ ~ 36K is on the edge of anisotropic scattering in higher 
order harmonics but also right in the middle of sin2(f) based anisotropic scattering. As we 
can see, the quantitative calculations for La — 214 compound contain a large uncertainty 
margin which is caused by a lack of a precise value of the specific heat jump at the phase 
transition and consequently of the electron density of states on the Fermi surface. The error 
in the experimentally estimated magnitude of AC is of the order of 36%. 

A similar analysis of Ni substituted samples shows that the pair-breaking effect in the 
unitary limit is stronger than observed experimentally. For the sake of comparison we 
present Ni doped F- 123 compound data0'01iiJli and our theoretical curves for the resonant 
scattering in Fig. 9. Less difference between the experimental and theoretical results is seen 
for La — 214.ii Although this comparison is suggestive for a weak potential scattering, the 
detailed calculations in the Born scattering limit giving the critical temperature dependence 
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on the residual resistivity^ are needed in order to draw more firm conclusions. 

Finally we discuss the electron irradiation experiments in F — 123.lll^ The results for 
low-energy (60-120 keV) incident electrons read from Fig. 12 of Ref. 21 are shown in Fig. 
10. The experimental points are in the range of the pair-breaking effect of the scattering 
potential with the anisotropy given by sm20. The initial suppression, however, is more 
gradual than the one obtained from the theoretical calculation and cannot be explained by 
the standard impurity pair-breaking mechanism. Unfortunately, we do not have any other 
set of electron irradiation data representing the change of the critical temperature with 
defect concentration. The analysis by Tolpygo et. al^ shows that the initial Tc suppression 
for low-energy (100 keV) electron irradiation reported by Legris et al.0 agrees with the one 
shown in Fig. 10. For higher-energy electronslll T^. is reduced initially about twice as much 
as in Ref. 21 which becomes in the range of the sin2(f) determined behavior for low defect 
concentration. 



VII. CONCLUSIONS 

We have studied the pair-breaking effect of the anisotropic impurity scattering in the 
t-matrix approximation for d^^-yi-wave and (c?a.2_y2 -|- s)-wave superconductors. The Born 
and the unitary limits have been discussed analytically and numerically. Although limited 
to the employed phenomenological model the conclusions we draw should be suggestive to 
a larger class of potentials capturing the feature of anisotropy. 

We have shown, that the effect of the anisotropic impurity scattering can be considered 
in four groups for the Born scattering and three groups in the case of unitary scattering 
which are determined by the form of the function / (k) defining symmetry of the scattering 
potential. In both scattering limits the strongest suppression of the critical temperature 
is caused by the impurity potential given by / (k) ~ cos20 which is in phase with the 
order parameter for the d^^-yi state and overlaps significantly with the {dx2_y2 + s)-wave 
state of a major d-wave component. The pair-breaking effect of this potential with a large 
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amount of anisotropic amplitude can be comparable with the one of the isotropic scattering 
in the Born limit and exceeds the s-wave impurity effect in the unitary limit. This issue 
is particularly important as it shows that a weak reduction of Tc due to anisotropy of the 
impurity potential proportional to the superconducting order parameter which follows from 
the weak-scattering modefl is not a general feature of the anisotropic impurity scattering. 
Another class is defined by / (k) ~ sin2(f). It leads to the lowest impurity pair-breaking 
effect in the case of the unitary scattering and second lowest for the Born limit. This kind 
of scattering is maximal in the direction of the nodes of the order parameter and it vanishes 
where the gap function has its maxima. Therefore, the effective scattering is minimized 
by the symmetry of the impurity potential. Any other function orthogonal to the (i-wave 
order parameter results in a rather universal T^. suppression and falls into the third group 
of the potentials. The pair-breaking in this case is less than that of the isotropic scattering 
but it exceeds the one of sm20 based potential, and in the unitary limit is very close to 
the isotropic scattering. Resonant scatterers in the p-wave channel lead to the same 
suppression as the (i-wave scattering given by sin2(f) function. In the Born scattering limit, 
however, the p-wave anisotropic scattering results in the lowest suppression. 

We have compared our results for the dx2_y2-wave superconductor with the experimental 
data assuming that the impurity potential is close to the unitary limit. Within the accuracy 
of our calculations the Zn atoms can be considered as the resonant scatterers in overdoped 
La— 214 and F— 123 compounds. The pair-breaking effect of the structure defects produced 
by the electron irradiation is also in the range of magnitude of the unitary scattering with 
the anisotropy of the impurity potential given by sm20 function. The theoretical results 
strongly depend on the value of the electronic density of states on the Fermi surface. Thus 
far, for some systems (like La — 214) there is only an estimation of this value with a large 
margin of uncertainty available. It opens a wide range of possible solutions and makes 
a proper analysis very difficult. Therefore, before any final conclusion about the impurity 
suppression of the critical temperature in the cuprates can be made, a much more accurate 
determination of the electronic density of states on the Fermi surface has to be done. 
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FIGURE CAPTIONS 



Fig. 1. Normalized critical temperature TjT^^ of the dj.2_j^2-wave superconductor as a func- 
tion of the normalized impurity scattering rate r'/(27rTc(,) in the Born limit for the impurity 
potential symmetry (solid curves from the top): sin24>, sin34> and higher-order harmonics, 
cos2(f). The amount of anisotropy in the potential is: (a) {1 — a) 0.17 {va/vi — 0.2), 
(b) (1-a) 0.33 {va/vi = 0.5), (c) (1 - a) = 0.5 {va/vi = 1.0), (d) {I - a) ^ 0.67 
{va/vi — 2.0), (e) {1 — a) ~ 0.83 {va/vi = 5.0) the higher-order harmonics split in this 
case into (from the top): sinA(l), cos3(l), sin3(l), cosA(l). The dashed line represents isotropic 
suppression of the critical temperature. 

Fig. 2. Shaded area represents the effect of the anisotropic impurity potentials given by 
cosl(f) and sinlcf) functions for 3 < Z < 10 and {1 — a) — 0.5 {va/vi = 1.0) in the dx2_y2-wave 
state in the Born limit. 

Fig. 3. The dx2_y2-weiYe state Tc suppression by the Born scattering p-wecve impurity po- 
tential with / (k) = sin4> (or cos4>) and the anisotropic scattering strength (1 — a) (va/vi) 
equal to (from the top) 0.83 (5.0), 0.67 (2.0), 0.50 (1.0), 0.33 (0.5), 0.17 (0.2). The isotropic 
scattering is shown with the dashed line. 

Fig. 4. Normalized critical temperature Tc/Tc^ for the {dx2-y2 -\- s)-wave superconductor 
((e) = 0.1) as a function of the normalized impurity scattering rate V in the Born 
limit for the impurity potential symmetry (solid curves from the top): cos0, sin0, sin24>, 
sin3(l) and higher-order harmonics, cos2(f). The amount of anisotropy in the potential is 
(1 — a) = 0.5 {va/vi = 1.0). The dashed curve represents the effect of the isotropic impurity 
scattering. 

Fig. 5. Normalized critical temperature T^/Tcg of the dx^^y^-wecve superconductor as a 
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function of the normalized impurity scattering rate V /{2t:Tcq) in the unitary hmit for the 
impurity potential symmetry (solid curves from the top): sin2(j) (and p-wave), sinicj) and 
higher-order harmonics, cos2(j). The dashed curve represents the effect of the isotropic im- 
purity scattering. 

Fig. 6. Normalized critical temperature T^/Tcq of the [dx^-y^ + s)-wave superconductor 
((e) = 0.1) as a function of the normalized impurity scattering rate V /{2'kTc^) in the uni- 
tary limit for the impurity potential symmetry (solid curves from the top): sin2(f> (and 
p-wave), sin34> and higher-order harmonics, cos2(j). The dashed curve represents the effect 
of the isotropic impurity scattering. 

Fig. 7. Critical temperature of Zn substituted optimally doped F— 123 samples vs impurity 
concentration per cent per planar Cu site, nzn- The data are taken from Refs. 13 (filled 
circles), 14 (open squares), 15 (diamonds), 16 (open circles), 17 (filled squares). The area 
between the curves of the same sort corresponds to the pair-breaking effect of the unitary 
scattering impurities of the following anisotropics: sm20 (solid), cos30 and higher-order 
(long dashed), cos20 (dashed). The isotropic impurity scattering is represented by the dot- 
dashed curves. 

Fig. 8. Critical temperature of Zn substituted La — 214 samples vs impurity concentration 
per cent per planar Cu site, nzn- The data are taken from Refs. 1 (filled squares), 3 (open 
squares). The theoretical curves in the unitary limit are drawn for two critical temperatures: 
(a) — 30ir, (b) Teg = ?>&K. The impurity potentials represented by different curves agree 
with the notation in Fig. 7. 

Fig. 9. Critical temperature of Ni substituted optimally doped 123 samples vs impurity 

concentration per cent per planar Cu site, riNi- The data are taken from Refs. 12 (open 
circles), 14 (filled squares), 15 (diamonds), 18 (open squares). The theoretically calculated 
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unitary scattering curves correspond to the potentials described as in Fig. 7. 

Fig. 10. Critical temperature of optimally doped F— 123 sample with in- plane oxygen va- 
cancies induced by 60 — 120 keV electron irradiation (Ref. 21) vs defect concentration per 
cent per planar Cu site, n(^o defects)- The theoretically calculated unitary scattering curves 
correspond to the potentials described as in Fig. 7. 
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FIG. 7. Fig.3: 
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FIG. 8. Fig.4: 
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FIG. 10. Fig.6: 
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FIG. 11. Fig.7: 
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FIG. 12. Fig.Sa: 
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FIG. 13. Fig.Sb: 
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FIG. 14. Fig.9: 
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FIG. 15. Fig.lO: 
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